The layered compound CaClFeP is an Arsenic-free high T c iron-pnictide 
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We first analyze why the iron pnictides with high T c superconductivity so far are As-based, by the 
Hund's rule correlation picture, then examine the P-based and Sb-based cases, respectively. Conse- 
quently, we propose that CaClFeP with ZrCuSiAs-type structure is an As-free high T c iron-pnictide. 
The subsequent density functional theory calculations show that the ground state of CaClFeP is of 
a collinearly antiferromagnetic order on Fe moments with structural distortion, resulting from the 
interplay between the strong nearest and next-nearest neighbor antiferromagnetic superexchange 
interactions bridged by P atoms, similar as the As-based pnictides. The other P-based pnictides 
are either nonmagnetic or magnetic but with weak exchange interactions. The Sb-based pnictides 
unlikely show high T c superconductivity because of the existence of robust ferromagnetic order. 

PACS numbers: 74.70. Xa, 74.20. Pq, 74.20.Mn 



The discovery of high transition temperature T c su- 
perconductivity in LaOFeAs by partial substitution of O 
with F atoms [l| stimulates the intense studies on the 
iron pnictides. Upon doping or high pressure, there are 
three types of iron pnictides reported to show supercon- 
ductivity, i.e. 1111-type ifeOFeAs (Re = rare earth) [l| 
or BaFFeAs and SrFFeAs 0, 122-type AFe 2 As 2 (A=Ba, 
Sr, or Ca) Q, and Ill-type BFeAs (B = alkali metal) 
All these pnictides share the same structural feature 
that there exist the robust tetrahedral layers where the 
Fe atoms are tetragonally coordinated by the As atoms 
and the superconduction pairing may happen. 

To our knowledge, all the iron pnictides showing high 
T c superconductivity so far are As-based. Here the high 
T c means close to or above the theoretical maximum 
value of 40K predicted from BCS theory [5j so that the 
corresponding superconductors are considered unconven- 
tional. Up to now, the highest T c reported in As-free 
pnictides is 17K shown in (Fe2P2)(Sr4Sc206) [6[. Is it 
possible to find a P- or Sb-based pnictide without As to 
show high T c superconductivity? Physically to answer 
this question essentially depends upon how we under- 
stand the mechanism underlying the high T c supercon- 
ductivity in the pnictides. Likewise, solving this question 
will greatly help us to understand the underlying mech- 
anism. Moreover, setting of this issue is very meaningful 
and important in practice because Arsenic strong toxic- 
ity sets the strict safety conditions to synthesize As-based 
pnictide samples, which confine the samples making in a 
few research groups worldwide. In this Letter, we show 
that CaClFeP with the ZrCuSiAs-type structure is an As- 
free high T c pnictide with the similar electronic structure 
and magnetic property as the As-based pnictides. 

It was universally found that all the As-based pnictides 
are in a collinear antiferromagnetic (AFM) order below 
a tetragonal-orthorhombic structural transition temper- 
ature 0, @] . Accordingly, the superconducting pairing in 
the As-based pnictides is now considered as being medi- 



ated by the spin fluctuations. Regarding the mechanism 
behind the structural and AFM transitions and their un- 
derlying relationship, there are basically two contradic- 
tive views, based on the itinerant electron picture |9j and 
the local moment picture [1CH12| respectively. From the 
very beginning, we proposed [12] the fluctuating Fe local 
moments with the As-bridged AFM superexchange inter- 
actions as the driving force upon the two transitions, ef- 
fectively described by the J1-J2 Heisenberg model. Here 
we would like to emphasize the twofold meanings in our 
proposal [Hj]: (1) there are localized magnetic moments 
around Fe ions and embedded in itinerant electrons in 
real space; (2) it is those bands far from rather than 
nearby the Fermi energy that determine the magnetic 
behavior of the As-based pnictides, namely the hybridiza- 
tion of Fe with the neighbor As atoms plays a substan- 
tial role. Here the formation of a local moment on Fe 
ion is mainly due to the strong Hund's rule coupling on 
the Fe 3d-orbitals, which is about 0.6-0.8 eV/Fe[l2j. In 
this sense, our proposal can be referred to as the Hund's 
rule correlation picture. We emphasize again that the 
Arsenic atoms play a substantial role in our physical pic- 
ture. Based on this Hund's rule correlation picture, we 
had successfully predicted that the ground state of a- 
FeTe is in a bi-collinear AFM order, which was confirmed 
by the later neutron scattering experiment 

The main factor competing with the Hund's rule cou- 
pling is the chemical bonding between the Fe 3d orbitals 
and As 4p orbitals. Meanwhile the hybridization between 
the Fe 3d orbitals and As Ap orbitals will bridge the su- 
perexchange interactions between the Fe magnetic mo- 
ments induced by the Hund's rule coupling. It is thus 
a dilemma to require strong superexchange interactions 
with robust Fe magnetic moments. Chemically, the spin- 
degenerated Fe 3d orbitals form stronger bonds with the 
As 4p orbitals than the spin-polarized Fe 3d orbitals do. 
For the As-based pnictides, inspection of the calculations 
[l2| shows that the spin degeneration pushes down the en- 
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FIG. 1: (Color online) CaClFeP with the ZrCuSiAs-type 
structure: (a) a tetragonal unit cell containing two formula 
units; (b) schematic top view of the FeP layer. The small 
red dashed square is an a x a unit cell, while the large black 
dashed square is a \[2a x \[2a unit cell. The Fe spins in the 
collinear antiferromagnetic order are shown by the red arrows. 



FIG. 2: (Color online) CaClFeP in the nonmagnetic state: (a) 
energy band structure, in which solid lines for the undoping 
while dotted lines for the 10%-replacement of Ca by La; (b) 
the Fermi surface; (c) top view on the Fermi surface, (d) and 
(e) correspond to the case of 10%-replacement of Ca by La. 



ergy levels raised by the Hund's rule coupling, and then 
strengthens the bonding between the Fe 3d orbitals and 
As 4p orbitals, which makes substantial energy gain to 
the nonmagnetic state. However, this is still not enough 
to completely overcome the Hund's rule coupling, in need 
of an additional energy gain of ~0.2 eV/Fe, as suggested 
by calculated total energy difference between the non- 
magnetic state and collinear AFM state [12]. Conse- 
quently, there are magnetic moments induced around Fe 
atoms with the AFM superexchange interactions bridged 
by the As atoms. 

In line of pnicogen atoms of P, As, and Sb, P has a 
much smaller atomic radius with a larger electronega- 
tivity than As while Sb has a much larger atomic ra- 
dius with a smaller electronegativity than As. Hence, it 
is very likely that in the most P-based pnictides there 
will form a very strong bonding between P and Fe so 
that the Hund's rule coupling is quenched. This is ver- 
ified by the calculations reported below that the most 
P-based pnictides are nonmagnetic. On the other hand, 
in the most Sb-based pnictides there will form a relatively 
weak bonding between Sb and Fe so that large moments 
are induced around Fe atoms with robust ferromagnetic 
spin fluctuations, which is also verified by the calcula- 
tions reported below. In contrast, Arsenic is in a proper 
balance to the dilemma. This explains why the iron pnic- 
tides with high T c superconductivity, discovered so far, 
are As-based rather than P-based or Sb-based. Now we 
deduce that if the in-plane lattice parameters a and b 
of a P-based (Sb-based) pnictide are more large (small), 
its electronic structure and magnetic property are more 
similar to the ones of the As-based pnictides since the 
P-Fe (Sb-Fe) bonding is more weakened (strengthened). 
Referring to the As-based pnictides, we further deduce 
that CaClFeP has the magnetic structure and property 
similar to the ones of the As-based pnictides, in great 
possibility to realize the high T c superconductivity. 

To verify the deduction, we carried out the density 
functional theory (DFT) electronic structure calcula- 
tions. In the calculations the plane wave basis method 



was used [15]. We adopted the generalized gradient ap- 
proximation (GGA) of Perdew-Burke-Ernzerhof [16| for 
the exchange-correlation potentials. The ultrasoft pseu- 
dopotentials [17j were used to model the electron-ion in- 
teractions. After the full convergence test, the kinetic 
energy cut-off and the charge density cut-off of the plane 
wave basis were chosen to be 600eV and 4800eV, respec- 
tively. The Gaussian broadening technique was used and 
a /c-mesh of 16 x 16 x 8 or 12 x 12 x 8 was adopted to 
sample the Brillouin zone of a x a x c or \/2a x y/2a x c 
supercell, respectively. In the calculations, the lattice 
parameters with the internal atomic coordinates were all 
determined by the energy minimization. CaClFeP is in a 
tetragonal layered structure with PA/nmm symmetry. A 
crystal unit cell consists of eight atoms with alternating 
FeP and CaCl layers along the c axis (Fig. QJa)). 

As analyzed above, our calculations show that Ill-type 
BFeP (B = Li or Na) and 122-type AFe 2 P 2 (A=Ba, Sr, 
or Ca) are all nonmagnetic semimetals since the strong 
Fe-P bonding quenches the Fe magnetic moments. Actu- 
ally, the lattice parameters a and b of these compounds 
are mainly determined by FeP layers, around 3.8 A. The 
large parameters a and b go to the 1111-type, in which 
the other layers strongly affect the lattice parameters. In- 
deed, we find the large lattice parameters (a=b=4.14 A, 
c=8.70 A ) for CaClFeP in magnetic phases. 

In Fig. [2] we plot the electronic band structure and 
Fermi surface of CaClFeP in the nonmagnetic state. 
There are three cylinder-like Fermi surface sheets, among 
which the two are of electron-type around M-A and the 
other one is of hole- type around T-Z. The volumes en- 
closed by these Fermi sheets give 0.19 electrons/cell and 
equally 0.19 holes/cell, i.e. 1.31 x 10 21 /cm 3 . The com- 
pound CaClFeP is thus a semimetal with a low car- 
rier density, between normal metals and semiconduc- 
tors. The density of states (DOS) at the Fermi en- 
ergy is 3.65 states/(eVx unit cell). The corresponding 
electronic specific heat coefficient 7 = 8.60 mJ/(K 2 * 
mol) and the Pauli paramagnetic susceptibility \v = 
1.49 x 10 -9 m 3 /mol. The 10%-replacement of Ca by 
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FIG. 3: (Color online) CaClFeP in the collinear-ordered anti- 
ferromagnetic state with the in-plane angle 7 = 90.4° (see Fig. 
1(a)): (a) energy band structure; (b) the Fermi surface; (c) 
top view on the Fermi surface. Notice that YX corresponds 
to the parallel- aligned moment line while FX' corresponds to 
the antiparallel-aligned moment line. 



La makes the electron-type and hole-type Fermi sheets 
more matched in shape, as shown in Fig. [2] (d) and (e). 

In order to explore the magnetic structures and clarify 
the exchange interactions between the Fe-Fe moments, 
we calculated three different magnetic states, the ferro- 
magnetic, checkerboard (Neel) AFM, and collinear AFM 
states, which respective energies are calculated to be 
(0.040, -0.044, -0.132) eV/Fe if the energy of the non- 
magnetic state is set to zero. The magnetic moment 
around each Fe atom is found about 2.0/i#. All these 
magnetic states can be self-consistently described by the 
following frustrated Heisenberg model with the nearest 
and next-nearest neighbor couplings J\ and J2, 



H 



Ji Si • Sj + J2 Si - s. 

(ij) 



(i) 



whereas (ij) and <C ij ^> denote the summation over 
the nearest and next-nearest neighbors, respectively. 
From the calculated energy data, we find that J\ ~ 
20.96 meV/S 2 and J 2 - 32.58 meV/S 2 (the detailed cal- 
culation is referred to Appendix of Ref. 12). The ground 
state of CaClFeP is thus in a collinear AFM order. 

In the collinear AFM phase, a small structural dis- 
tortion is found to further make energy gain by 1.2 
meV/Fe, which changes the angle between two prin- 
cipal axes in ab-plane 7 from 90° to 90.4° (see Fig. 
[TJb)), equivalently, the lattice unit cell slightly expands 
along the spin-antiparallel direction and slightly shrinks 
along the spin-parallel direction. Here we emphasize that 
our calculations show that the driving force upon this 
tetragonal-orthorhombic structural distortion is nothing 
but the magnetic interaction. More specifically it is the 
superexchange interaction J 2 that induces the Fe-spin 
collinear AFM order and drives the structural distortion 
to break the rotational symmetry, as found in LaOFeAs 
[12|. Moreover, the anti-parallel alignment for the Fe 
moments between the neighboring FeP-layers was found 
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FIG. 4: (Color online) Calculated total and projected density 
of states at the five Ye-3d orbitals around one of the four Fe 
atoms in the collinear antiferromagnetic state unit cell. 



TABLE I: Calculated exchange couplings J1-J2 (Fig. D^b)) 
with measured superconductivity transition temperature T c . 





Ji (meV/Fe) 


J 2 (meV/Fe) 


Tc (K) 


LaOFeP 


0.52 


7.83 


7[6] 


(Fe 2 P 2 )(Sr 4 Sc 2 06) 


3.74 


12.24 


17[21] 


BaFFeP 


0.62 


11.20 


? 


CaClFeP 


20.96 


32.58 


? 



almost degenerate with the parallel alignment. 

From Fig. [3j we notice that the Fermi surface is more 
distributed along ZR r or YX' than along ZR or YX di- 
rection in Brillouin zone for CaClFeP in the collinear 
AFM state, similar as in LaOFeAs and BaFe 2 As 2 [l2|,[l8[. 
This suggests that the spin-antiparallel direction is more 
conductive than the spin-parallel direction in real space. 
This ant i- intuitive finding was confirmed by the recent 
experiment on the detwinned BaFe 2 As2 [Hj]. 

Fig. H] projects the density of states onto the five 3d or- 
bitals of Fe in the collinear AFM of CaClFeP. We see that 
the five up-spin orbitals are almost filled while the five 
down-spin orbitals are nearly uniformly partially filled. 
This shows that the P atom-imposed crystal field split- 
ting is small. The formation of Fe magnetic moments is 
thus due to the Hund's rule coupling, which is a universal 
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feature found in all the As-based pnictides 

Other 1111-type P-based pnictides were also calcu- 
lated. They all show the similar electronic structures 
and magnetic properties but with much smaller exchange 
couplings, as reported in Table HI These compounds have 
smaller lattice parameters (about 4 A) than CaClFeP. 

We have known that all the As-based pnictides have 
the similar electronic structures and magnetic properties. 
Upon doping or high pressure, the AFM long-range or- 
der is suppressed. The remanent AFM fluctuations are 
now considered to be responsible for the superconducting 
pairing in pnictides [12j, similar as in cuprate supercon- 
ductors. We had systematically calculated the exchange 
couplings J 1 and J 2 for the different As-based pnictides. 
And we found that phenomenologically the superconduc- 
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FIG. 5: (Color online) Calculated magnetic energy versus the 
angle between the two sublattice spin orientations, from 
the ferromagnetic order (0 = 0°) to the Neel AFM order (6 — 
180°). (a) NaFeAs and (b) NaFeSb. 

tivity transition temperature T c is proportional to J 2 [HI • 
Theoretically, the two-band t-J\-Ji model also indicates 
that the T c in pnictides is proportional to the J2 . 

As we have just shown, CaClFeP takes the similar elec- 
tronic structure and magnetic property as the As-based 
pnictides and the exchange couplings J\ and J2 of it are 
as large as the ones of BaFe2As2 [Hj]. We thus pre- 
dict that the high T c (likely over 40K) superconductivity 
can be realized in CaClFeP. Meanwhile, for LaOFeP and 
(Fe2P2)(Sr4Sc20e) the smallness of calculated J\ and J2 
is consistent with the observed low T c (Table UJ). 

We have performed the DFT calculations on Sb-based 
pnictides, including Ill-type, 122-type, and 1111-type 
compounds. We find that all these Sb-based pnictides 
are similar with each other in electronic structures and 
magnetic properties. In magnetic states, they show rich 
magnetic orders like the checkerboard AFM, bi-collinear 
AFM, and collinear AFM orders, however, there are 
much larger moments formed at Fe atoms (~ 3/ie/Fe) 
in the Sb-based pnictides than in the As-based pnictides 
because of the strong Hund's rule effect and the weak 
Fe-Sb bonding, as analyzed above. Besides, even though 
the collinear AFM order is the most stable energetically, 
there exists a stable ferromagnetic order in the Sb-based 
pnictides, independent of the AFM fluctuations. This is 
essentially different from the As-based pnictides, in which 
a ferromagnetic order is never stable. 

Fig. [5^a) shows that for NaFeAs the ferromagnetic 
order is a maximum energetically, and the magnetic en- 
ergy versus the angle between the two sublattice spin 
orientations is well described by the nearest neighbor ex- 
change Heisenberg model. The calculations also show 
that CaClFeP behaves likewise. In contrast, for NaFeSb 
the ferromagnetic state is a local minimum in the energy 
landscape, lower than the nonmagnetic state (Fig. [5jb)). 
And the variation of the magnetic energy with the angle 
9 severely deviates from the exchange Heisenberg model. 
This shows that the ferromagnetic state has a different 
exchange interaction origin than the AFM states, which 
is likely related to the possible orbital ordering. This 
ferromagnetic state can be stable until applying a high 



pressure of 25 GPa. Since such ferromagnetic correla- 
tions are against the superconducting pairing induced by 
the AFM fluctuations, we don't think there is the high T c 
superconductivity realized in Sb-based pnictides. To our 
knowledge, no Sb-based pnictides are reported to show 
high T c superconductivity in experiment. 

In summary, by the Hund's rule correlation picture we 
may well understand the iron pnictides, in which the dif- 
ferent magnetic structures and properties result mainly 
from the competition between the Hund's rule effect on 
Fe 3d-orbitals and the Fe-P, Fe-As, or Fe-Sb bonding. In 
search for an As- free high T c pnictide, we show that Ca- 
ClFeP is an AFM semimetal with the strong P-bridged 
nearest and next-nearest neighbor AFM superexchange 
interactions, which give rise to the collinear AFM order 
on Fe moments in the ground state, similar as LaOFeAs 
and BaFe2As2- Upon doping or high pressure, CaClFeP 
is predicted to show high T c superconductivity. The other 
1111-type P-based iron-pnictides are also AFM semimet- 
als but with weak exchange interactions, while both Ill- 
type and 122-type P-based iron-pnicitdes are nonmag- 
netic semimetals. We also show that the Sb-based pnic- 
tides unlikely show high T c superconductivity. 
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